Exact solution for a non-Markovian dissipative quantum dynamics 
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We provide the exact analytic solution of the stochastic Schrodinger equation describing an har- 
monic oscillator interacting with a non-Markovian and dissipative environment. This result repre- 
sents an arrival point in the study of non-Markovian dynamics via stochastic differential equations. 
It is also one of the few exactly solvable models, for infinite dimensional systems. We compute the 
Green's function; in the case of a free particle, and with an exponentially correlated noise, we discuss 
the evolution of Gaussian wave functions. 
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Non-Markovian quantum dynamics are rapidly becom- 
ing a popular field of research. From the theoretical point 
of view, there is a increasing interest in understanding 
the behavior of physical systems beyond the Markov ap- 
proximation [IH1]. This interest is triggered by novel 
experiments [SJ E] , which do not find explanation within 
the standard Markovian description. Moreover, people 
are already speaking of possible future applications, for 
example in describing energy transfer in biological sys- 
tems 0. 

There are different approaches to the description of 
(Markovian and non-Markovian) open quantum dynam- 
ics [MTTj. among which stochastic differential equations 
(SDEs). In this approach, one modifies the Schrodinger 
equation by introducing new stochastic terms, which 
mimic the effect of an external agent (typically, the en- 
vironment). This method has particular advantages: it 
allows for a description in terms of state vectors, instead 
of density matrices; it helps in finding exact analytic so- 
lutions for the dynamics, as in the case presented here; 
in some cases, it allows for faster numerical simulations. 

A particularly relevant example of SDE is the follow- 
ing [IMS]: 



d 
dt 



V~Xqwt — Xq 2 



(1) 



here, q is the position operator of the particle, H its quan- 
tum Hamiltonian, A is a positive coupling constant, and 
Wt is a Gaussian white noise. The equation does not pre- 
serve the norm of the state vector; this is not a problem, 
as at any time the vector can be appropriately normal- 
ized. This equation plays a key role in several fields. 
Within decoherence theory, it represents one of the pos- 
sible unravelling of the Joos-Zeh master equation [20] : 
within collapse models, it represents the simplest descrip- 
tion of a quantum particle undergoing spontaneous col- 
lapses in space [12] : within the continued measurement 
theory, it describes a system whose position is continu- 
ously measured Q]. In the following, we will focus our 
attention on collapse models, though our results apply 
also to the other fields. 

If one tries to identify the noise w t with a physical field, 



one has to cope with the fact that it presents two unphys- 
ical features: it is an infinite-temperature field, due to the 
lack of dissipative terms, and has a white noise spectrum, 
i.e. each frequency contributes with the same weight. 
These features can be disturbing, and people worked in 
improving the model. The dissipative generalization of 
Eq. 0) has been developed in [21]; in the new dynamics, 
the noise is such that each system thermalizes to a fixed 
temperature, that of the field. The Markovian behav- 
ior has been generalized by studying the non-Markovian 
version of Eq. , where the white noise is replaced by a 
general Gaussian (colored) noise [HI [22l [23] . 

In this Letter, we present a novel result [24]: the gen- 
eralization of Eq. both to finite-temperature (dissipa- 
tive) and colored (non-Markovian) noises, and its solu- 
tion. The equation is the following: 



dt 



-2VXq 



t c 

dsD(t,s) — 



(2) 



The new terms depending on the momentum operator p 
account for dissipation, whose strength is determined by 
the positive constant //. The temperature of the noise 
is determined by the relation: T = H 2 (4mK B fj,) . The 
integral term of Eq. 0, which involves the whole past 
history of the system, accounts for the non-Markovian 
behavior of the dynamics: D(t,s) is the time correla- 
tion function of the noise, and r— denotes a functional 
derivative with respect to w s . 

Eq. plays a central role, because it reproduces all 
previous models, under the appropriate limits. When 
D(t,s) —> S(t — s), it reduces to the Markovian dissipa- 
tive model analyzed in 21 . When fj, — > 0, it reproduces 
the non-dissipative and non-Markovian model presented 
in [TTJ [53]; when both limits are taken, it reduces to 
Eq. 0. 

We succeeded in finding the general solution of Eq. , 
for an harmonic oscillator [24] . which we will describe 
here below. Taking into account that only for equations 
containing terms which are most quadratic in p and q, 
a general method for finding the solutions is known, our 
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equation represents the most general equation that can 
be analytically solved (modulo physically unimportant 
terms). This matches with the fact that, from the phys- 
ical point of view, the model represents an arrival point, 
since it involves both a non-Markovian and dissipative 
noise. Note that the model is similar to that of Hu, 
Paz, Zhang, describing non-Markovian quantum Brow- 
nian motion [8]. 

Solution. The solution of Eq. (pi is given in terms of its 



Green's function G(x, t; xq, 0). In [35] it has been shown 
how G(x,t;xo,0) can be written via the path integral 
formalism as follows: 

G(x,t;x ,0) = X V[q]e s ^, (3) 

Jq(0)=x 

where S[q] is a non standard Action. For our model, it 
takes the expression: 



S[q} = 



q^-m\^q s q s -^Q 2 q^+mV\fj,q s - A s q s +q s / dr B(r, s)q r - 2mXfi q s / drD(s,r)q r 
* Jo Jo 



with Q = uj — A /i , where u> is the proper frequency of the oscillator, and 



, (4) 

(5) 
(6) 

An important issue here arises. Since S[q] is time-non-local, the standard Lagrange formalism cannot be applied; 
therefore we resorted to a time-non-local variational formalism in |26| . The path integration in Eq. Q was performed 
using the mid-point formulation of the Feynman polygonal approach |27| . The calculation is long, and is duly reported 
in |24j : here we present only the final solution, which is an exact result. The Green's function reads: 



A s = ihV\w s + rriX 3/2 fi 2 w s + 2mA 3/ V 2 / drD(r, s)w r , 

Jo 

B(r, s) = (2mA V + 2ihX) D(r, s) + AmX 2 ^ 2 [ dr'D(r, r')D(s, r') , 



G{x,t;x Ol 0) = 



exp 



-A t xl — A t x 2 + B t x x + C t x + V t x + £ t 



(7) 



At = k[ / t (0) — Xfx — 2X/J, / dsD(0,s)f(s) 



Bt = k[ f t (t) - g t (0) + 2A/i / ds D{0, s)g(s) 



2iithtu(t) 

where u(t) is a suitably defined function, unimportant for the subsequent analysis, At = —k (gt(t) — A/i), k — Wr, and 

(8) 
(9) 
(10) 

(11) 



Ct = -k(ht(0)+2VXfiw + VXu f w s f t (s)ds+[ — f t (s)ds - 2X(i f dsD(0,s)h(s) 
V -'o Jo m Jo 

V t — k ( h t (t) + 2y/Xfi w(t) — V X^i [ w s g t (s)ds — [ — g t (s)ds) , 
V ./o Jo m J 

t 



-k ^/Xn J w s ht(s)ds + J 



A. 



ht(s)ds — Xfi / Wg ds ] , 



where the differentiation is always done with respect to the variable between parenthesis. Defining the following 
integro-differential operator: 



I[e(s)] ■= j e(s) + — (n 2 + AX^D(s, s)) e(s) - ^ 



dr B(r, s)e(r) 



1 



dr B(s,r)e(r) — mXfi I dr 



dD{r, s) 
dr 



e(r) — mX/i / dr 



dD(r, s) 
~ds 



e(r) . 



(12) 



I 

one can prove that /(s), g(s), and h(s) solve the following with suitable boundary conditions 
intcgro-diffcrential equations: 



The Green's function ^ has some important features. 
I[f(s)} = mXnD(0,s) I[g(s)} = 0, (13) The first one is that it has a Gaussian structure. This 

implies, first of all, that Gaussian states evolve into Gaus- 

(14) 

sian states. It implies also that the density matrix has a 



m^/j A(s) 
= — - — w(s) — 
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Gaussian structure, because it is the average (an integral) 
with respect to a Gaussian noise, of states which evolve 
under a Gaussian-preserving dynamics. A second prop- 
erty is that the coefficients A, A and B are deterministic, 
while only C, V, and E depend on the noise Wt and thus 
evolve stochastically. The fact that the first three coef- 
ficients do not depend on the noise wt guarantees that 
the spread of a wave function evolves deterministically 
in time. From the point of view of collapse models, this 
is an important property, because it guarantees that any 
wave function collapses in space. Finally, one can check 
that, as expected, in the limit ^ — > the solution of the 
non-Markovian model is reproduced [llj . while in the 
limit D(t, s) —> 6(t — s) one recovers the Markovian dis- 
sipative model [21]. Considering both limits one obtains 
the Green's function of Eq. (JlJ. It is worthwhile stress- 
ing that the proof of all these properties relies on the 
fact that Eqs. (|7|)-(|l~4|) represent the analytical solution 
of Eq. j2). 

The explicit expression of the Green's function depends 
on the solution of the three integro-differential equa- 
tions ( 13 1-( 14 ). These equations can be exactly solved 



only for some correlation functions D(t,s). In order to 
give a quantitative example of the features of the model, 



we have solved Eqs. ( 13 1-( 14 ) for the exponential correla- 



tion function and have studied the behavior of the spread 
of Gaussian wave functions. 

Exponential correlation function. Wave function for a 
free particle. Let us consider an exponential correlation 
function 



D(t,s) = ( 7 /2)e-^- s \ 



(15) 



which, besides allowing for an explicit solution of the 
equations, also represents a good physical example of a 
correlation function with a finite correlation time, given 
by 7 _1 . Using Eq. (15), one can show that the integro- 



differential equations (13)-(14) can be transformed into 

One can solve 



fourth-order differential equations [28 
these equations for a generic quadratic potential (in par- 
ticular, the harmonic oscillator [23]). For simplicity, here 
we consider a free particle, whose equation of f(s) reads: 

J(s) - ( 7 2 + AV - 2A /n ) f(s) 

2ih\ 1 2 s 



+ (4AVV + 
The general solution is 



f(s) = 0- 



(16) 



ft(s) = }Jft,k sinhw fc s + g t ^ k cosh v fc s] , (17) 



k=l 



where the coefficients ft,k, 9t,k are determined by the 
boundary conditions, and v\, vi are the two roots of the 



characteristic polynomial associated to Eq. ( 16 1 



ui, 2 = v / ((7-A/i) 2 ±C)/2, 
C= -1/(7 - A/i) 4 - 16A V7 2 



8ihXj 2 



(18) 
(19) 



Two of the four boundary conditions are ft (0) = 1 and 
ftit) = 1, while the other two can be determined using a 
standard procedure explained in [23l [28] : 



/(0) = AV/(0) + A//7 / die 
Jo 



-7« 



/(/) 



+ A/17 + 



ihX-f 



III 



(lit 



-7* 



/(/) 



(20) 



/(f) = AV/(*) - AA7 / dle-^-»f(l) 



3A 2 /i 2 7 i/LVy\ 



■J. 



dle-^-^fH) 



AV7 



di e - 7(t+J) /(0 



(21) 



These four boundary conditions allow to determine the 
exact form of the coefficients ft^k and gt,k] their expres- 
sions are too long to be written here, but can be easily 
worked out e.g. with Mathematica® . The same proce- 
dure here described can be applied also in order to find 
the analytic expressions for g(s) and h(s). All the details 
of the calculation are reported in [23] ■ 

Having found the explicit expression for every coeffi- 
cient of the Green's function ^ , one can analyze the time 
evolution of wave functions. Particularly interesting, and 
easy to analyze, are Gaussian states, whose form — as pre- 
viously discussed — is preserved by the dynamics. Ac- 
cordingly, a wave function of the type: 



4> t {x) = cxp[-a t x 2 + (3 t x + 74] 



(22) 



is solution of Eq. ([2]), and the coefficients at, fit and 74 
evolve in time as follows: 



at = At 



It = 7o + £t + 



B 2 t 

4(a +A t ) 

(Ct + M 2 



(3t = 



Ct+Po 
4(a + A t ) 



4(a +A t ) 



(23) 



As anticipated, one can easily see that the evolution of 
the spread both in position and in momentum is deter- 
ministic since the parameters A, A and B do not depend 
on the noise. In particular, we focus our attention on the 
behavior of the spread in position: cr(f) = l/2y/af (the 
apex r denotes the real part) in the case of the exponen- 



tial correlation function (15). 



Fig. [T] shows the comparison of the evolution of the 
spread according to the four versions of the model 
(white [2H], white dissipative [5T], non- white [TT], non- 
white dissipative) . First of all one can see that the be- 
havior is qualitatively the same for every model: the 
spread decreases in time, reaching an asymptotic finite 
value. From the quantitative point of view one can see 
that both in the non-Markovian (dot-dashed line) and 
in the dissipative (dashed line) models the shrinking of 
the wave function is slower than in the white noise case 
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FIG. 1. Time evolution of the spread in position a(t) of a 
Gaussian wave function, in the case of an exponential corre- 
lation function for the noise. The initial spread is cr(0) = 1 
m. The other parameters have been chosen as follows: m = 1 
kg, A = 3 x 10 24 m" 2 s" 1 , fi = 1CT 30 m 2 (dissipative models), 
7 = 10 -4 s _1 (non-white models). Time is measured in sec- 
onds, distances in meters. For very small values of <r(0) the 
wave function spreads out instead of shrinking, reaching the 
asymptotic value displayed in the figure (see also [317]). 



(solid line). As a consequence, when these two effects 
are combined (dotted line), the spread decreases even 
more slowly. This is easy to understand. In the Marko- 
vian case, all frequencies contribute to the collapse of 
the wave function, while in the non-Markovian case the 
high-frequency components are suppressed. This slows 
the collapse process. Secondly, a finite temperature noise 
is less energetic than an infinite temperature noise, there- 
fore in the dissipative case the collapse is slower than in 
the non-dissipative case. 

From the point of view of collapse models, an inter- 
esting physical question is whether a field with 'typical' 
cosmological values for the temperature (T ~ 2.73K, 
that of the CMBR) and for the spectrum (cutoff at 
7 - 10 10 - 10 n Hz, that of the CMBR, the relic neutrino 
background, and the relic gravitational background), can 
collapse the wave function efficiently. This question has 
been answered positively in [31] : a noise with typical 
cosmological features can collapse the wave function fast 
enough to guarantee the emergence of classical properties 
at the macroscopic level. 
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